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Abstract 

The fluctuations and correlations of matrix elements of cross sections are investigated in open systems that are 
chaotic in the classical limit. The form of the correlation functions is discussed within a statistical analysis and 
tested in calculations for a damped quantum kicked rotator. We briefly comment on the modifications expected 
for systems with slowly decaying correlations, a typical feature in mixed phase spaces. 
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1. Introduction 

The photo-dissociation cross section as deter- 
mined by Fermi's golden rule contains a combi- 
nation of final density of states and transitional 
matrix elements. The statistical properties of this 
cross section in a situation of a chaotic dynamics 
will thus be determined by the statistical proper- 
ties of both the density of states and the matrix 
elements. Moreover, if there is a classical under- 
lying dynamics, both will be connected, in a suit- 
able semiclassical limit, to the properties of the 
classical dynamics. Our aim here is to present a 
few numerical and theoretical considerations con- 
nected to these observations. 

We focus on the form of the correlation func- 
tion of photo- absorption cross sections. Using 
a random matrix theory approach Alhassid and 
Fyodorov [1] found a correlation function that 
consisted of a Lorentzian and a derivative of a 
Lorentzian with respect to the width. The first 
term is familiar from the analysis of Ericsson fluc- 
tuations in nuclear physics [2,3] and the second 
one from the correlation function of the Wigner 
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time delay [4] . The form of the correlation func- 
tion can now be made plausible if within a sta- 
tistical model for the cross section the fluctua- 
tions in the density of states and in the transition 
matrix elements are independent, as discussed in 
section 2. 

A semiclassical analysis of this correlation func- 
tion [5] shows that the relative weight of the two 
contributions depends on the ratio of the fluctu- 
ations of the observable to the mean. In case of 
a single initial state, random matrix theory fixes 
this ratio to universal numbers. However, in situ- 
ations an incoherent superposition of initial states 
contributes to the cross section, variations in the 
relative weight are possible. This is illustrated 
within a numerical analysis for a damped kicked 
rotator in section 3. 

The semiclassical connection also suggests cer- 
tain modifications in the correlation functions if 
the classical decay is not purely exponential[6]. In 
particular, in situations with mixed phase space 
an algebraic decay is expected[7]. The modifica- 
tions in the correlation functions include a slower 
decay and the formation of a cusp at the origin, 
depending on the exponent of the decay law. In 
section 4 we propose a model for the form factor 
and analyze some of the consequences. 

We conclude with a brief summary in section 5. 
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2. Matrix element correlations within ran- 
dom matrix theory 

The quantity we focus on is the density of states 
weighted by the matrix elements of the observ- 
able, 



p A {E) = Y j A ll 8 v {E-E ll ), 



ft 



(1) 



where the sum runs over the eigenstates of 
the system having eigenvalues E^. The photo- 
absorption cross section is proportional to this 
expression if the observable A contains the pro- 
jection onto the initial state and the dipole oper- 
ator. The expectation value of A can be written 



as A^ = A 



<L4 M , where A is the mean and SA^ 



is the random fluctuation around its mean. Sta- 
tistically we assume 

(5A n ) = 0, and (5A n 5A m ) = a A 5 n , m . (2) 

The function 5^ (x) is a Lorentzian function with 
half width parameter r], normalized so that as 
77 > it approaches a Dirac-5. We take the 
same value of rj for all the eigenstates implying 
a uniform damping. The mean density of states 
is simply 



p A = (pa(E)) e 



dE ,^ NA 
-^-Pa(E) = — , (3) 

where N is the number of levels and B is the 
energy width of the subset of the spectrum, over 
which the average (■ ■ -)e is calculated. 

The normalized autocorrelation function of the 
fluctuations of the density of states, Spa(E) = 
Pa(E) — pa is defined as 

(S PA (E + e)Sp A (E)) E 



C(e)-. 



Pa 



(4) 



(na) 



dE_ 
~B 



PA (E + e) PA (E) - 1. 



Inserting the definition (1) in (4) we obtain 

C{£) = (ll) 2 ^ ^ + A[5A » + 5Av) 

+5A li 5A v ]g 1 ,(e,E lt ,E v )-l, (5) 
where we have inserted the shorthand notation 
9ri {e,E^E u ) = (5 ri {E-E li + e)5 r ,(E-E v )) E .(6) 



With the assumption that matrix elements and 
resonances are uncorrelated [8], averaging over 
the SA's eliminates two terms, 



C(e) 



Y,9v(e,E^E,)-l (7) 



N 



( B 



2 



+ [na) J2( sa ^ a M £ ,e„e u ). 

In the second term we will utilize (2). Also we will 
split the first double sum in diagonal (p = v) and 
non-diagonal (p ^ v) parts. For sufficiently small 
r\ the diagonal term is the autocorrelation func- 
tion of Lorentzian that also yields a Lorentzian. 



C(s) = (1 + 

B 

N 



+ TT 



Za 

A 2 
2 



El 

N 

9v( £ > E^E,,). 



(5 v (E-e)5 v (E)) E -l 
(8) 



After performing the averaging in the first term 
and in the sum together with the definition in (6), 
one arrives at 



C{e) = A^(e)-1 



(9) 



where A = B/N. In the sum one can recog- 
nize the appearance of the two-level correlation 
function in the limit of 77 — > 0. For strongly over- 
lapping resonances, i.e. when 77 3> A the above 
expression reduces to 



C{e) oc 



V 



+ ■ 



1 



A 2 V 2 + e 2 2n(r] 2 +e 



2\2 



(10) 



Thus the correlation function is characterized by 
two terms, a Lorentzian and a derivative of a 
Lorentzian with respect to the broadening pa- 
rameter 77. The weight of the first term comes 
from the fluctuations of the observable A. In the 
case when A = is a projection on the basis 
state This quantity in random matrix the- 
ory (RMT) [9] is (J3 + 2)/(3, where j3 = 1,2, and 
4 for the different universality classes the system 
belongs to (orthogonal, unitary, and symplectic, 
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respectively). Hence we recover the correlation 
function derived in Ref. [1] for channels with uni- 
form resonance width. 

However, if A vanishes, the normalization by 
A is not possible and the second term, which 
comes from the correlation function of the den- 
sity of states, cf. (5) and (6), disappears and the 
correlation function becomes a pure Lorentzian. 
These findings are thus in accordance with what 
has been argued on the basis of semiclassical pe- 
riodic orbit theory [5]. 

3. Correlations in the quantum kicked ro- 
tator 

To illustrate the above calculations we calcu- 
late the correlation function of the matrix element 
weighted density of states for a damped kicked ro- 
tator. 

We consider the statistical properties of observ- 
ables that are projections onto a subset of the 
basis states 

a= E HH ( n ) 

where I(m) is a subset of size m < N of the ba- 
sis set. The fluctuations of the matrix elements 
(p,\A\p,) over the eigenstates /U of the system de- 
scribe the cross section fluctuations of the exci- 
tation of the system from an initial state I(m) 
to the final state \i [1,10,11,5]. A possible dipole 
operator has been absorbed into the definition of 

Im>- 

The model system we considered is a quantum 
kicked rotator with a kicking potential 

V(4>) = fc(cos</>- sin 20). (12) 

It is known [3,6] that this model belongs to the 
unitary universality class since the second term in 
the potential breaks the conjugation symmetry 

We have diagonalized the unitary one-step evo- 
lution operator U at a value of the classical kick- 
ing strength, K, where complete ergodicity was 
expected (K = 7 [12]). The size of the system was 
fixed to N = 201. The matrix element weighted 
density of states was defined as 



M0 = 5>|iM 



M-r/2 



)~\ (13) 



The mean value of this observable is simply A = 
m/N where 1 < m < N is the number of states 
contributing in the projection. We have calcu- 
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Figure 1. The correlation function C(x) as a 
function of the eigenvalue separation in units of 
mean level spacing. The data points are obtained 
for observables of projections over m basis states 
indicated in the legend. The continuous curves 
are fits of the form (15) fitting a and 7 indepen- 
dently. The width of the resonances was uniform 
throughout the spectrum: T = A. 

lated the correlation function (4) applied for our 

case 



C(x) = (Sp A {<P + x)5p A ((j))) 4n 



(14) 



i.e. averaging (. . is done over the eigenphase 
spectrum extending over B — 2ir. The variables 
4> and x are measured in units of the mean level 
spacing A = B/N therefore pa = 2nA/N. The 
arguments of [5] and the RMT arguments of the 
previous section shows, that the correlation func- 
tion is composed of a Lorentzian and a derivative 
of a Lorentzian 



C(x) 



7 



1 



2 2 
7 Z — x z 



~f 2 +x 2 27r( 7 2 + a; 2 ) s 



(15) 



In (15) a = a A /(A 2 T H ), with the Heisenberg 
time Th = N. The width 7 should come out 
to be the damping 77 of (10), in appropriate units. 

This expression can be compared to numerical 
results on the kicked rotator. In Fig. 1 we plot the 
correlation function obtained for observables that 
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are projections extending over m — 1, 3, 10, 50 ba- 
sis states. The continuous curves in the figure are 
fits of the form (15) allowing the two parameters 
a and 7 to vary. The classical estimate for a is 
[5,6] 



A 



a = 



A 2 N 



1 

— = oc — 

AN m 



(16) 



The variance a\ follows its classical value for low 
values of p = A = m/N [6]. We obtained a fitted 
value of 7 = 1.11 r (instead of T) independent of 
manda= 0.284, 0.135, 0.042, and 0.009 for m = 
1, 3, 10, and 50, respectively. The dependence on 
m is close to the one predicted by (16) at least 
for m > 1. 

4. Open systems and slowly decaying cor- 
relations 

The correlations described in the previous sec- 
tions can be derived also using the following phe- 
nomcnological procedure [1]. The two-level corre- 
lation function of a closed chaotic system, Co(e), 
consists of a Dirac-<5 at the origin and a smooth 
function decreasing to zero for large level separa- 
tions. This is the Fourier transform of the form 
factor K (t). From now on for sake of simplicity 
we restrict ourselves to the unitary universality 
class, i.e. we write the RMT form factor in its 
standard form [9], K (r) = 1 — 6(r), with 



&(r) = (l-|r|)0(l-|r|) 



(17) 



where 9(x) is the step-function. Time r is mea- 
sured in units of the Hciscnbcrg time T H , there- 
fore energy separation e, in units of mean level 
spacing A. According to RMT the correlation 
function is 

/oc 
drKoir) cos(27rer) 
-OO 



= 6(e) 



sin(7re) 



7T£ 



(18) 



Uniform damping, i.e. opening up the system, can 
be introduced by multiplying the form factor with 
an exponential decay. 



where T measures the relaxation rate in units of 
the Hciscnbcrg time. In fact T = T c /Th, where 
T c is the relaxation time. When T < 1 (r > 1) 
the relaxation happens over a time scale longer 
(shorter) than Tjj, hence produces correlations 
over energy scales lower than (beyond) mean level 
spacing. The resulting form factor for different 
values of T = 0.0, 0.1, 1.0 and 5.0 are plotted in 
Fig. 2. The Fourier transform of (19) results in 



C{e)=C l {e) + C 2 {e) 
where 

/oc- 
dre- r|r| cos(2tt£t), 
-OO 



(20) 



(21) 



C2(e) 



|r|)e" r|r| cos(27T£t). (22) 



The first term C\(e) is a Lorentzian from the 
broadening of the (5-function and the second term, 
the convolution of the Lorentzian with the two- 
level cluster function, leads, in the limit of large 
r, to the derivative of the Lorentzian. The cor- 
responding correlation functions C (e) are plotted 
in Fig. 3. 




(19) 



Figure 2. The form factor for different uniform 
damping constants T = 0.0, 0.1, 1.0, and 5.0. 

These considerations led us to extend the above 
approach in another direction. Let us assume 
that the quantum return probability, the form 
factor, is damped slower than exponential, i.e. 
in an algebraic fashion. The correlation function 
in this case will contain signatures of the slowing 
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Figure 3. The correlation function, C(e) (see 
Eq. (20)), of an open system for different uniform 
damping rates V = 0.0, 0.1, 1.0, and 5.0. The 
inset shows C^e) (22). 

down of the classical dynamics, a behavior that is 
expected to be prominent in systems with mixed 
phase space [13]. 

Such an algebraic damping may be introduced 
in the phenomenological ansatz 



K(t) = K (t)(1 + ctY 



(23) 



where c = Th/T c is the ratio of the decorrela- 
tion time T c compared to the Heisenberg time 
Th ■ As in the case of exponential damping, for 
c < 1 (c > 1) the slow decorrelation of classi- 
cal trajectories due to the presence of the hier- 
archy of stable islands [13,14] occur on a time 
scale that is longer (shorter) than the Heisenberg 
time that produces correlations over energy scales 
lower than (beyond) mean level spacing. By fix- 
ing one of the parameters a and c the variation of 
the other results in changes in different parts of 
the correlation function. Similarly as in (20) we 
find 

C{e) = C 1 {e) + C 2 {s) (24) 
where 

J-oo (1 + W 

C 2 (s) = -j\ T j^L corner). (26) 

We have plotted the correlation function for the 
case when a = 4/3. This type of correlations can 
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Figure 4. The form factor with algebraic damping 
(a = 4/3) for different ratios of the Heisenberg 
time Tjj compared to the decorrelation time T c , 
c = T H /T C = 0.0, 0.1, 1.0, and 5.0. 

be expected to occur in systems with mixed phase 
space [14]. In Fig. 5 the curves demonstrate that 
a low value of c results in only a slight modifi- 
cation of the correlations while large c, i.e. when 
T c < Th the algebraic decorrelation results in 
changes both in the e < 1 and e > 1 regimes. 
This is in contrast to previous expectations [14] 
that the dynamics occuring over time scales up 
to Th should show up in the correlation function 
on energy scales beyond mean level spacing only. 

To show the similarities and differences be- 
tween correlation functions obtained for exponen- 
tial and power law damping we plot the value 
at e = 0, i.e. the variance of the fluctuating ob- 
servable under investigation. In Fig. 6 we can 
see that for T c > Th the two types of damping 
functions yield the same correlation time depen- 
dence, oc (Th/Tc)- 1 . In the case of T c < Th the 
exponential damping produces a (Th/T c )~ 2 de- 
pendence and the power law version a (Th/Tc)- 11 
(in the case shown a = 4/3). 

An even more striking difference between the 
fast and slow decorrelations shows up in the low- 
£ behavior of the correlation function. As pointed 
out already by Lai et al. [15] non- hyperbolic sys- 
tems produce a cusp in the correlation function 



C(e) ~ C(0) - C a e a ~ 
where C a = baC-^- 1 



(27) 

is a positive constant and a 



is the exponent of the classical return probability 
P(t) ~ t~ a . This behavior is to be contrasted 
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Figure 5. The correlation function, C(e), of a sys- 
tem containing algebraic decorrelation with expo- 
nent a = 4/3 (see Eq. (24). The curves stand for 
different ratios of the Hcisenberg time Th com- 
pared to the decorrelation time T c , c = Th/T Ci 
with c = 0.0, 0.1, 1.0, and 5.0. The inset shows 
C 2 (e) (26). 



with the hyperbolic case when the low-£ behavior 
is expected to be 

C(e) ~C(0)-C 2 e 2 , (28) 

where C2 ~ r~ 2 is a constant depending on the 
exponential relaxation rate. In Fig. 7 we show 
the cusp at low-£ in the case of power law decor- 
relation and the parabolic behavior for exponen- 
tial damping. The curves represent especially the 
cases when Th /T c > 1, i.e. when the classical cor- 
relations appear at time scales below the Heisen- 
berg time. In this case one expects deviations at 
£ > 1. The cusp at low-e, however, is present for 
any value of c. 

We would like to emphasize that the apparent 
contradiction between the behavior of the corre- 
lation function for power law decorrelation (27) 
and that of rcfs. [13,14] resides probably in the 
difference of the classical return probability func- 
tions considered. We took the modification of 
the RMT form factor (23) while in rcfs. [13,14] 
the classical return probability of P(t) ~ t z was 
used. 



Figure 6. The value of C(e) at e — for exponen- 
tial (continuous) and power law damping (dashed 
curve). The ratio Th/T c is T for exponential and 
c for power law damping. The power law curve 
was obtained for the exponent of a = 4/3. 
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Figure 7. The correlation function for low-£ for 
exponential (continuous) and power law damp- 
ing (dashed curve) for two different values of 
the decorrelation time T c /Th = 1, and 5. The 
power law curve was obtained for the exponent 
of a = 4/3. 

5. Final remarks 

The variations in the correlation functions for 
the different situations are not very large and per- 
haps difficult to detect. This seems to apply in 
particular to cases with an algebraic decay, espe- 
cially in view of the fact that the exponents are 
not universal and might be clouded by a distri- 
bution of algebraic decay laws. The modifications 
due to variations in matrix elements and thus a 
stronger emphasis of the derivative part of the 
correlation function could perhaps be achieved in 



incoherent superpositions of cross sections from 
different initial states. Analysis of experimental 
data in this direction seems worthwhile. 
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